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Abstract 



o 

m . 

We point out that the world sheet swept by the confining string in 
presence of dynamical quarks can belong to two different phases, de- 
pending on the number of charge species and the quark masses. When 
it lies in the normal phase (as opposed to the tearing one) the string 
breaking is invisible in the Wilson loop, while is manifest in operators 
O i[ composed of disjoint sources, as observed in many numerical experi- 

ments. We work out an explicit formula for the correlator of Polyakov 
loops at finite temperature, which is then compared with recent lattice 
data, both in the quenched case and in presence of dynamical quarks. 
The analysis in the quenched case shows that the free bosonic string 
model describes accurately the data for distances larger than ^0.75 
fm. In the unquenched case we derive predictions on the dependence 
of the static potential on the temperature which are compatible with 
the lattice data. 



e-mail:gliozzi, provero@to.infn.it 



1 Introduction 



The confining interaction of a pair of static sources in non Abelian gauge 
theories is mediated by a thin flux tube, or string, joining the two sources. 
At large distance, the energy of the flux tube is proportional to its length, 
hence the static potential grows linearly. 

When matter in the fundamental representation is added to the system, 
one expects that this potential flattens at some distance r , where the string 
breaks to form pairs of matter particles which screen the confining potential. 
The broken string state describes a bound state of a static color source (the 
fixed end of the string) and a dynamical matter field (the free end of the 
string). A similar screening is expected even in pure Yang- Mills theory 
for adjoint sources, since in this case the adjoint string may terminate on 
dynamical gluons. 

Despite many efforts, these string breaking effects have proved elusive 
in standard analyses of large Wilson loops, both in QCD with dynamical 
fermions [Q] and in pure Yang-Mills theory with adjoint sources in 3+1 
dimensions [|j and in 2+1 dimensions ||. 

On the contrary, clear signals of string breaking have been observed 
in studies where the basis of operators has been enlarged in order to find 
a better overlap to the true ground state, following a method originally 
advocated in Ref. In this way, fundamental string breaking has been 
found in SU(2) Higgs model in 2+1 || and 3+1 @ dimensions, and adjoint 
string breaking in 2+1 SU(2) pure gauge theory |7], ||[. 

The outcome of these analyses is twofold. On one hand the Wilson loop 
appears to have in general very poor overlap on the broken string state, 
the only exception being observed up to now in the 2+1 dimensional SU(2) 
gauge theory with two flavors of staggered fermions ||] . On the other hand 
it turns out that the operators with a good overlap with the broken string 
state have as a common feature the presence of two disjoint source lines. 
This is not only true in the above mentioned cases at zero temperature, 
but it is also evident in the recent observation of string breaking at finite 
temperature QCD with dynamical fermions, where the static potential is 
extracted from the correlator of two disjoint Polyakov loops jlO| . 

In this work we will show that in the usual string description of confine- 
ment [|ll]] one can find a simple explanation of such a relationship between 
the overlap properties and the number of disjoint static sources. In partic- 
ular we shall derive an asymptotic functional form of the unquenched finite 
temperature potential which is then successfully compared with the QCD 
data of Ref. @. 
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A solvable prototype of string breaking |L^ ] suggests that the world sheet 
swept by the string in its time evolution can exist in two different phases, 
known as tearing and normal phases. The question of the existence of a 
finite overlap of the Wilson operator with the broken string state can be 
reformulated (see section [2]) as the problem of finding the phase of the string 
world sheet. Indeed in the former phase the world sheet is torn by large holes 
corresponding to pair creation of charged matter, hence the Wilson loop is 
expected to fulfil the perimeter law at large distances. On the contrary, 
in the normal phase the string breaking is balanced by the inverse process 
of string soldering, so that adding dynamical matter to the system yields 
simply a renormalization of the string tension. In particular the Wilson 
loop behaves exactly like in the quenched model and no macroscopic string 
breaking is visible. 

However, a simple topological argument shows that the situation dras- 
tically changes when one considers, instead of the Wilson loop, operators 
constructed out of disjoint static sources. Besides the usual contribution to 
the confining potential there is a new term which is necessarily absent in 
the quenched case and which survives at large source separations. This term 
makes string breaking effects easily detected on the lattice. 

The remaining sections are organized as follows. In section || we re- 
view some known bosonic string formulae about the asymptotic infrared 
behaviour of gauge operators and derive the new term contributing to cor- 
relator of two Polyakov loops in the unquenched case. Since there are no 
recent studies on the string effects in the quenched approximation of QCD, 
we devote section |] to an analysis of the quenched data of Ref. JlPf] . It 
turns out that the bosonic string accurately describes the lattice data for 
distances larger than ~ 0.75 fm. In section || the analysis is extended to the 
unquenched case where we check that string predictions on the temperature 
dependence of the potential are compatible with the lattice data of Ref. [1C]. 
Finally in section ^ we draw some concluding remarks. 



2 A topological argument 

Let us start by considering the standard Wilson loop. In absence of charged 
dynamical fields this loop acts as a fixed boundary of the surface associated 
to the world sheet swept by the confining string. When matter fields are 
added, any number of holes of any size may appear on this surface, reflect- 
ing the pair creation of dynamical charged particles. If these are light, the 
holes behave like free boundaries. The total string contribution to the Wil- 
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son operator is then the sum over all possible insertions of multiconnected 
boundaries of the world sheet. It can be written diagrammatically as a loop 
expansion (see Fig. 1). 

There are two quantities controlling the number and the size of these 
holes. One is the mass of the matter fields (it is easier to create lighter 
particles). The other is the number of charge species. For a theory with Nf 
flavors and N c colors there is a factor of NfN c for each hole. Likewise in 
the adjoint string every hole is accompanied by a factor of N% — 1. 

For light particles and NfN c large enough, we expect that the string 
world sheet is dominated by configurations with a large number of relatively 
small holes compared with the size of the Wilson loop. 

These holes do not influence the area law of the Wilson loop and their 
effect can be absorbed into the renormalization of the string tension, ac- 
cording to a string mechanism originally proposed in Ref. Under the 
circumstances there is no way to observe the string breaking through the 
study of large Wilson loops. This may explain why this phenomenon has 
been so elusive both for fundamental and adjoint string in any dimension. 
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Figure 1. Loop expansion of the rectangle. 



However it is worth noting, as mentioned in the introduction, that there is 
a solvable matrix model [jl^] describing idealized random surfaces with dy- 
namical holes (or open strings embedded in zero dimension), where besides 
the above described "normal" phase dominated by holes of small size there 
is also another phase, characterized by a spontaneous tearing of the surface, 
due to the formation of large holes growing with the size of the surface (see 
Fig. 2); as a consequence we expect that the associated Wilson loop should 
fulfil the perimeter law.[] 

What is the phase of the world sheet of the confining string? The ob- 
served extremely poor overlap of the Wilson loop with the broken string 
state seems to suggest that it belongs to the normal phase^], with a possible 

1 There is also a third phase separating the above two, where both the world sheet (the 
glue) and the holes (the charged matter) are large and in competition. This phase could 
play a role at the deconfining temperature. 

2 Of course the numerical data cannot exclude the possibility that such a poor overlap 
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Figure 2. Configuration of a torn surface 



exception for the SU(2) gauge theory in 2+1 dimensions with dynamical 
fermions |J. Perhaps in this case the quarks are not light enough. 

For the operators with more than one connected source line the loop 
expansion has a topologically different form. Take for instance a pair of 
Polyakov loops of a gauge theory with dynamical fermions at finite temper- 
ature. From the point of view of the effective string, they are the two fixed 
boundaries of a cylinder (see Fig. 3). The loop expansion splits into two 
different series, because starting at two-loop level there are configurations 
in which the two fixed boundaries are not connected by the string world 
sheet. 

If the system is in the normal phase the sum over the loop insertions 
can be still absorbed into the renormalization of the string tension, but now 
the loop expansion is the sum of two different terms. That with only one 
connected world sheet has no overlap to the broken string state but decays 
exponentially with the distance between the fixed boundaries, thus only the 
term with two disjoint pieces survives, which of course is expected to have 
a large overlap to the broken string state. In the case of light matter fields 
this term can be written as a square [Zdn(R , L t )] 2 , where Zdn denotes the 
contribution of a cylindric world sheet with fixed (or Dirichlet) boundary 
conditions (b.c.) on the Polyakov line of length Lt and free (or Neumann) 
b.c. on the other side (representing a dynamical charge) placed at a mean 
distance R Q from the Polyakov loop. Clearly R Q fixes the scale of the string 
breaking. Owing to the simple geometry of this term, we shall explicitly 

is a finite size effect and that larger Wilson loops will eventually reveal a spontaneous 
tearing of the world sheet. 
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Figure 3. Loop expansion for the cylinder. 



evaluate its functional form in the next section. 

Similar considerations apply to operators used to see string breaking 
at zero temperature. In such a case the fixed boundary is not necessarily 
closed because it may terminate on charged fields. The rule to associate a 
string world sheet to these operators is simply to connect the end points 
of the fixed boundaries (Wilson lines) with paths associated to the free end 
of the string, representing the world lines of the dynamical particles. For 
instance, using the diagrammatic representation of Ref. ||, the map to the 
string picture of other operators used to determine the static potential in 
the SU(2) Higgs model is drawn in Fig. 4. 
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Figure 4. Map to string world sheet of operators used to determine the 
static potential in the SU(2) Higgs model. The filled circles correspond to 
Higgs sources, the thick lines represent Wilson lines or fixed boundaries, 
the thin lines free boundaries. 
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3 String formulae 



The main role of the string picture of confinement is to fix the functional 
form of gauge operators in the infrared limit. In particular it is known jlj] 
that, according to the bosonic string model, a rectangular R x L Wilson 
loop should behave asymptotically as the partition function of d — 2 free 
two-dimensional bosonic fields describing the transverse oscillations of the 
world sheet with fixed b. c: 



(W(R,L)) o, e —RL-p(R+L) 



7?(t) 



d-2 

2 



where d is the spacetime dimension, a the string tension, p the perimeter 
term and 77 is the Dedekind function 

00 

V(r) = q 1/24 J] C 1 " O > 1 = exp(2^r) . (2) 

71=1 

This in turn gives the static potential 

V{R) = - ^lim I \og(W(R, L)) =oR- + . . . (3) 

where dots indicate subleading terms. These predictions have been tested 
to high accuracy for the Z2 gauge model in d = 3 in |l5| . Good agreement 
with the string model predictions was also found for SU(N) gauge theories 



in d = 3 |U and for SU(3) ind = 4 |7|. 



Similar considerations apply to the confining phase of finite temperature 
gauge theories. Here the relevant observable is the correlator of Polyakov 
loops. Since the lattice is finite and periodic in the time direction, the 
Polyakov loop correlation will be described by the partition function of d — 2 
free bosons on a cylinder with periodic b.c. in the time direction and fixed 
b.c. on the two loops. The partition function is therefore: 

—aoRLt 

(P(0)P + (R)) oc Z DD (R, L t ) oc ( 4 ) 

where o"o is the zero temperature string tension at the same coupling, Lt is 
the lattice extension in the time direction, that is the inverse temperature, 
and now r = ijfe. Sec. || will be devoted to a comparison of Eq. @) with 
lattice quenched QCD data at finite temperature. 
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We have seen in Section ^ that in presence of dynamical matter the 
Polyakov correlator contains an additional term [Zdn(R , Lt)] 2 , where Zjjn 
is the partition function on a cylinder with fixed b.c. on the Polyakov 
loop side and free b.c. on the other side, corresponding to a dynamically 
generated loop. The functional form of Z^n is similar to that of Eq. (B), 
but now the the r/ function is replaced by the function 

OO . j 

V(T ) = q~ 1/48 11(1 - qM) , r = y~, q = exp(2vrir ). (5) 

n=l ° 

The form af the two functions r] and fj can be simply understood as the 
inverse of the partition function of the normal modes of vibration of the 
string. When both ends are fixed the allowed frequencies are uj n = -^n and 
these generate the infinite product of Eq. (Q); if one end is free we have 
instead £j n = ^(n — |) which generate the infinite product of Eq. (||). The 

prefactor q^i in Eq. (|2|) is directly related to the zero point energy E of the 
string 

h-K ^ hm Hit 

n=l n 

where ^ denotes the ^-function regularized sum. It is known that within 
such a regularization one can safely handle the series as they were finite 
sums, then we can write [18] 

/ / / / 1 

5> = ^2n + £(2n-l) = 2C(-l)+2^>--) ; (7) 

hence 



2' 

n=l n=l n=l n=l 



^2 2.R ^ v 2 7 48i? v 7 

n=l n=l 

This gives in turn the prefactor of rj. The identity r){r) = ^fy^y allows us 
to write 

Z DN (R ,L t ) = e-°° R ° L ^ nW 



lv(r /2)\ ■ (9) 

Thus, assuming that the world sheet is in the normal phase, we would expect 
that the Polyakov loop correlation function in presence of dynamical quarks 
has the following effective string description: 

{P(0)P + (R)) = c 1 Z DD (R ,L t )+c 2 [Z DN (R ,L t )\ 2 



e ~aoRLt 



v( T o) 
v(t /2) 



2(d-2) 

(10) 
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where N c and Nf are respectively the number of colors and flavors of light 
quarks. The constants c\ and C2 and the length scale R Q are not predicted 
by the model and will have to be determined numerically. Notice that the 
second term in Eq. ( |Tq| ) is just a constant at any fixed temperature: therefore 
the actual predictions of the model concern the temperature dependence of 
the potential. In Sec. |5| we will compare these predictions with lattice QCD 
data. 



4 The static potential in quenched QCD at finite 
temperature 

If no dynamical quarks are present, the effective string prediction for the 
Polyakov loop correlation is given by Eq. (||). For the static potential we 
have in d = 4 

V{R) = -i-log(P(0)P+( J R)) = aoR+^-logr, {^) (11) 

We will be interested in the 2R > Lt region, in which Eq. (|Tl|) is conveniently 
rewritten in the equivalent form 

=°° R -^ t + h^ 2 i + iX l < 1 - e " 4 '" s/i ') (12) 

1 n=l 

Therefore at fixed temperature 1/Lt and asymptotically for large R we have 
a temperature dependent string tension 

<y{L t ) = cr - ^ ( 13 ) 

and a logarithmic term, plus exponentially suppressed subleading terms. It 
is perhaps worth stressing that there is no 1 / R term in this regime []. 

We want to compare the prediction Eq. (|l^) with Monte Carlo data for 
quenched QCD at finite temperature. First let us stress some important 
points: 

3 The Liischer term ir/12R appears in the opposite limit 2R -C Lt, as can be seen by 
rewriting Eq. ([H]) in the other equivalent form 

n— 1 
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1. The free string picture we are using is an effective infrared description: 
we expect it to describe the long distance behavior of the potential 
(for a detailed discussion of this point in the zero temperature case 
see Ref. [§1|). Therefore we expect to find a minimum distance i? s tring 
above which Eq. (^) describes the data accurately. Obviously there 
will be a distance i?ii ne ar above which the data are well described also 
by the simple linear behavior V ~ aR. The string model will be 
confirmed if i? s tring is sensibly smaller than -Rii nea r- 

2. The string model gives the finite temperature potential in terms of 
the zero temperature string tension <jq. Therefore the result of the 
fit of the finite temperature potential should be compared to the zero 
temperature string tension to confirm the string picture. 

3. The free string picture must break down at temperatures near the 
deconfinement point where string interactions are believed to become 
important (see e.g. Ref. jOJ). Therefore we need to use lattice data 
at temperatures not too close to the critical one. 

4. The main effect of string interactions is to renormalize the value of 
the string tension do while leaving the functional form Eq. (^) ap- 
proximately unchanged. Therefore we expect string interactions to 
introduce a systematic error in the evaluation of the zero temperature 
string tension from finite T data. 



Very precise data were obtained by the authors of Ref. |T(J at three 
temperatures, with T/T c ~ 0.8, 0.88 and 0.94. The first of these data 
samples is ideal for our purpose. We fitted those data to Eq. (|ll]) and, for 
comparison, to a simple linear behavior V ~ aR. Notice that these are 
both two parameter fits: the string contribution to the potential does not 
contain any adjustable parameters. Let us define -Rstring as the value of R 
above which the reduced x 2 of the fit is < 1, and -Rimear as the corresponding 
distance for the linear fit. The data are taken at (3 = 3.95 and Lt = 4. We 
obtain, in terms of the lattice spacing a, 

-Rstring = 3.3a (14) 

with x 2 = 0.74, and 

^linear = 4.5a (15) 

If we tried to fit the data for R > i? s t r ing to V = aR we would obtain 
X 2 = 8.6. 
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The fit with the string model potential and R > -Rstring gives 

ao = 0.22a" 2 (16) 

This should be compared with the string tension at zero temperature at the 
same value of the coupling (5 = 3.95. Extrapolating the data of Ref. [jHJ 
to our P we find ctq = 0.24(l)a~ 2 : we conclude that the systematic error 
discussed above is still rather small at T/T c = 0.8 (the statistical error on 
<7o given by our fit is two order of magnitude smaller and hence not very 
meaningful) . 

It is interesting to express i2 s tring in physical units: using 

a = (420 MeV) 2 = 4.41 fin -2 (17) 

we obtain 

^string ~ 0.75 fm (18) 
while using the string tension as the length scale 

ffoJ&ring ~ 2.4 (19) 

If one studies the data samples closer to the critical temperature, at 
T/T c ~ 0.88 and 0.94, one clearly sees the effect of string interactions dis- 
cussed above: while the fit of the potential with the free string model po- 
tential is always very good, the systematic error in the determination of 
o"o increases: At T/T c ~ 0.88 the fit to Eq. ([T^) gives o"o ~ 0.18a" 2 when 
the Wilson loop value is @ 0.197(8)a~ 2 , and for T/T c ~ 0.94 we obtain 
a ~ 0.13a" 2 instead of 0.173(5)a" 2 . 



5 String breaking 

In this section we compare the prediction Eq. ([To| ) to lattice data taken 
by the authors of Ref. |l(| in finite temperature QCD with two flavors of 
staggered dynamical quarks. For the static potential Eq. ([!(]) predicts, for 
d = 4, 



V(R) = - -log(P(0)P+( J R)) 

Lit 



-aoRLt 



T lo g Cl 



log 



7](t) 2 
aoRLt 



+ (N c N f f c 2 e- 2aoRoLt 



v(to/2) 



1 4 



7](t) 



(N c N f ) 2 c(R ,L t ) 



+ A 



(20) 
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where A = — log c\ and 



c(R ,L t ) = ^l e -^oRoL t 
ci 



v(ro/2) 



(21) 



The predictive content of the model lies in the dependence of c on Lf. 
To verify this prediction we fit the Monte Carlo data for the static potential 
to the expression 



V(R) 



los 



- (to RLt 



7](t) 



(N c N f ) 2 c 



A 



(22) 



and verify that the dependence of c on L t is actually described by Eq. (|2ip. 
We used data taken by the authors of Ref. [[HJ at L t = 4 and four values of 
(3 ranging from (3 = 5.1 to (3 = 5.28. The results of the fits of the four data 



samples to Eq. (21) are reported in Tab. 1, in lattice units: 






a a 2 


(N c N f fc 


x 2 


5.10 


0.594(80) 


0.0151(68) 


1.58 


5.20 


0.475(27) 


0.0297(47) 


1.11 


5.25 


0.407(19) 


0.0653(72) 


1.02 


5.28 


0.308(21) 


0.189(25) 


0.21 



Table 1: Fit of the lattice static potential to Eq. 
P considered 



for the four values of 



To verify Eq. (|2l|) we need to trade the (3 dependence at fixed Lt for a 
temperature dependence, by determining the /3-dependent lattice spacing 
a((3). We chose to use the values of do derived from our fit to evaluate 
a{(3) (remember that <To in E q. JlC| ) is the zero temperature string tension). 
Therefore fixing do using Eq. ( |l7l) we obtain the lattice spacing at the four (3 
values and therefore the inverse temperature in fm . Using [2C] T c /y/aQ = 
0.436(8) we can then determine the ratio T/T c for each (3. These values 
are reported in Tab. 2. With the exception of the highest temperature, for 
which we cannot trust the free string picture for the reasons mentioned in 
the previous section, the values of T/T c are in good agreement with the ones 
quoted in Ref. JiO|]. 

Now we can compare our results with Eq. (^l]). The latter contains two 
free parameters, namely the length scale R Q and the ratio c<ijc\. It is clearly 
safer to disregard the highest temperature point; unfortunately in this way 
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p 


L t (fm) 


T/T c 


5.10 


1.468(96) 


0.744(52) 


5.20 


1.313(28) 


0.832(28) 


5.25 


1.215(28) 


0.899(27) 


5.28 


1.057(36) 


1.033(39) 



Table 2: Inverse temperature in fm and the ratio T/T c as computed from the 
values of ao for the four beta values. 



we have to perform a two parameter fit of three data. Therefore the best we 
can hope at present is to show that the Monte Carlo data are compatible 
with our model, while more data at temperatures not too close to T c would 
be needed for a more conclusive test of the model 



The fit of the three values of c with Eq. (21) gives a satisfactory reduced 
X 2 of 0.64. The length scale R is 

R = 0.71(16) fm (23) 

The fact that this length scale turns out of the same order of magnitude as 
the typical string breaking scale is of course encouraging. The constant ^ 
is affected by an error which is bigger than the value of the constant itself: 
a = 4.2(6.4). 

If we push the model beyond its natural limits and use the parameters 
given by our fit to predict the value of c at the highest temperature, we 
obtain for f3 = 5.28 the value (NfN c ) 2 c = 0.20, again in good agreement 
with the numerical result. Even if this last result must be taken with all 
the necessary caution, it is nevertheless a significant hint that the model is 
realistic. 

We can conclude that the free string model of string breaking Eq. (|i~0| ) 
is compatible with lattice result for finite temperature QCD, and certainly 
reproduces at least the qualitative features of the phenomenon. More data 
at moderate temperatures are certainly needed to draw a definite conclusion 
about the correctness of the model. 



6 Concluding remarks 

In this paper we conjectured that the world sheet swept by the confining 
string in presence of matter fields belongs to the normal phase, characterized 
by a large number of microscopic holes produced by the matter fields, whose 
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net effect is just to renormalize the string tension. This must be contrasted 
with the tearing phase, in which large holes dominate and drive the string 
tension between static sources to zero. 

This conjecture has two main consequences, both amenable to numerical 
verification: 

• String breaking cannot be detected by studying the large distance 
behavior of Wilson loops: even in presence of dynamical matter fields 
the Wilson loop expectation value should behave at large distance like 
in the quenched case, that is according to Eq. (|l|). 

• On the contrary, when one considers operators made of disjoint static 
sources a simple topological argument shows that even if the world 
sheet is in the normal phase string breaking can be easily detected. 

Both of these are in agreement with most of the rather large body of lat- 
tice data on the subject that has been published recently O-[10|. The only 



exception is (2+l)-dimensional SU{2) theory with two flavors of staggered 
fermions, where string breaking effects in the Wilson loop have been ob- 
served ||. We suggest that the mass of the quarks used in Ref. |J is large 
enough to make the world sheet belong to the tearing phase. Indeed, the 
very same theory does not show any sign of string breaking in the adjoint 
Wilson loop 0, U where the "matter" is massless by definition. 

Based on this conjecture, we have proposed an effective string model 
of string breaking for the simplest example of an operator constructed out 
of disjoint sources: the correlator of Polyakov loops in finite temperature 
QCD. While the predictions of this model appear to be compatible with the 
lattice data, more data would be needed for a conclusive assessment of its 
validity. 

One can envisage a number of numerical experiments to further check the 
above conjecture and the effective string model. Perhaps the simplest could 
be the analysis of the Polyakov correlator in the adjoint representation at 
finite temperature. In this case one can work directly in the pure Yang Mills 
model and test accurately Eq. (|io|) as we did for Eq. (Q) in the quenched 
case. 

Also, it would be interesting to verify that the Wilson loop in presence 
of dynamical quarks actually behaves like in the quenched case, that is 
according to Eq. ([!]). It is however important to note that fuzzed Wilson 
loops should be avoided in this kind of check, because the string quantum 
fluctuations produce a strong shape dependence which is out of control in 
fuzzed loops. The r] function of Eq. (pi) accounts for these fluctuations only 
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in the case of a rectangular shape. Likewise, analyses based only on the 
asymptotic form of the potential given in Eq. (^) should be taken with 
caution, because the logarithmic correction generated by the yR term is of 
the same order of magnitude as the Liischer term within the lattice size of 
current simulations. 
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